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The microscopic dynamics and aging of a soft thcrmoscnsitivc suspension was investigated by look- 
ing at the thermal fluctuations of tracers in the suspension. Below and above the glass transition, 
the dense microgel particles suspension was found to develop an heterogeneous dynamics, featured 
by a non Gaussian Probability Distribution Function (PDF) of the probes' displacements, with an 
exponential tail. We show that non Gaussian shapes are a characteristic of the ensemble-averaged 
PDF, while local PDF remain Gaussian. This shows that the scenario behind the non Gaussian 
van Hove functions is a spatially heterogeneous dynamics, characterized by a spatial distribu- 
tion of locally homogeneous dynamical environments through the sample, on the considered time 
scales. We characterize these statistical distributions of dynamical environments, in the liquid, 
supercooled, and glass states, and show that it can explain the observed exponential tail of the 
van Hove functions observed in the concentrated states. The intensity of spatial heterogeneities 
was found to amplify with increasing volume fraction. In the aging regime, it tends to increase as 
the glass gets more arrested. 



1 Introduction 

Dense colloidal suspensions are a powerful model system to study the glass and jamming tran- 
sitions. With increasing volume fraction, the suspension dynamics may experience a slowing 
down of several orders of magnitude, as well as a change of qualitative behaviour, from vis- 
cous to clastic, at reasonable experimental time scales, with no significant structural changes in 
the suspension— This colloidal glass transition is analogous to the glass transition in liquids 
and polymer melts and to the jamming transition in granular materials^. Because the colloidal 
particles are larger than molecules, these systems provide slower accessible times scales than in 
atomic and molecular glasses, and the possibility to be probed with optical techniques including 
microscopy and dynamic light scattering—. Hard spheres suspensions have been extensively in- 
vestigated for their relative simplicity in modelling the glass transition—. Recently, soft colloidal 
suspensions have also attracted attention. They can be flexibly designed with various particles 
softness and interactions - attractive or repulsive - allowing for a wider variety of behaviour in 
mimicking the glass transition in traditional glasses^. One of the most investigated soft systems 
consists of aqueous suspensions of thcrmoscnsitivc microgels, made with amphiphilic cross-linked 
poly(N-isopropylacrylamid) polymer (pNIPAm)— . These particles can be designed with various 
softnesses, depending on the cross-linkers density—. Moreover and not least, the particles radius 
in the suspension, and therefore their volume fraction, can be reversibly tuned with temperature, 
which provides a unique way to explore the phase behaviour on the same sample. This flexibil- 
ity has been used to study the phase and behavioural transitions of colloidal suspensions^, such 
as crystal formation— ^ or meltin g 11 ! 12 , organisation in constrained geometry—, and glass and 
jamming transitions^ - — . 
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Colloidal glasses are intrinsically out-of-equilibrium systems which age : their physical proper- 
ties depend on the time elapsed since they were formed. When approaching the glass transition, 
colloidal suspensions exhibit heterogeneous dynamics^ - — , also observed in other glass forming 
systems^— and in simulations^— . It is characterized by non Gaussian tails of the van Hove 
correlation functions, and is now admitted to be a characteristics of the glass and jamming tran- 
sitions-^^. On a theoretical point of view, many interpretations were proposed to capture the 
heterogeneous dynamic s 25 i 34 i 35 , Experimentally, a connection with cage motion rearrangements 
with string like motions of high mobility, was established by different groups^—. Recently, it 
was related to a double Gaussian in a soft thermosensitive glass, suggesting two distinct dynamic 
populations of particles^. In another recent investigation in thermosensitive suspensions, mul- 
tiple regimes of dynamic heterogeneities were found with increasing volume fraction—. While 
widely investigated, the scenario of dynamic heterogeneity in glass forming materials, still remains 
elusive. 

In this paper, we report an experimental investigation of the heterogeneous dynamics in a ther- 
mosensitive soft spheres suspension when approaching the glass transition from below (increasing 
volume fraction), and beyond. Our system consists of a suspension of pNIPAm microgel parti- 
cles, with soft repulsive potential interactions, with temperature as the external control parameter 
of the volume fraction. We first investigate the aging dynamics of the soft glass by looking at 
the thermal fluctuations of the tracers. We then characterize the dynamic heterogeneity in the 
suspension by computing the ensemble-averaged and the local Probability Distribution Functions 
(PDF) of the probes' displacements. We characterize the statistical distributions of the dynamical 
environments explored by the tracers. At equilibrium in the supercooled regime, and in the glass 
state, we show that it can explain the exponential tail of the van Hove functions. Finally, the in- 
tensity of the spatially heterogeneous dynamics was investigated with increasing volume fraction, 
below and above the glass transition. 

2 Material and Methods 

2.1 Soft microgel particle suspensions 

Our system consists of a suspension of micron size cross-linked pNIPAm microgel particles, syn- 
thesized as follow : 20 mg of the cross-linker BIS (N,N'-Methylenebisacrylamide, Polysciences, 
Inc.), 20 g of NIPA (N-isopropylacrylamide, Polysciences, Inc.), and 400 ml of deionized water 
are loaded into a special three-neck flask equipped with a stirrer, thermometer, and a gas inlet. 
The resultant mixture is stirred, heated at 82 °C, and bubbled with dry nitrogen for 10 min to 
remove dissolved oxygen. A solution of 200 mg of APS (Ammonium Persulfate) dissolved in 1 ml 
of deionized water is then added to the mixture to start the polymerization reaction. The mixture 
is continuously stirred at 82 °C for 30 min and then allowed to cool down to room temperature. 
The resultant particles are centrifuged and resuspended in deionized water few times to remove 
unreacted monomer, homopolymers, and other salts. 

The microgel particles arc swollen at low temperatures, T < 30 °C, and collapsed above 
T > 35 °C, as shown in Figure [T] Below T = 32 °C, where we investigate the system behaviour, 
they have soft repulsive interactions. The hydrodynamic radius of the particles is a monotonously 
decreasing function of temperature, the LCST being around T = 32 °C. In colloidal hard spheres 
suspensions, the phase behaviour is controlled by the volume fraction. Concentrated colloidal hard 
spheres suspensions are known to enter an out-of-equilibrium glass state above a glass transition 
volume fraction around $ g ~ 0.58i Here, the volume fraction can be monitored with the temper- 
ature, allowing for an exploration of the phase diagram on the same sample. The polydispersity, 
of the order of 10% for our samples, as well as the particles softness, suppresses the crystallisation 
at high volume fractions (obtained at low bath temperatures). Starting from a liquid state and 
increasing volume fraction, the suspension was found to enter a glass state, characterized by aging 
over long time scales. 

A concentrated stock suspension of pNIPAm particles was prepared. Polystyrene beads with 
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a diameter of 1 /im were dispersed in the pNIPAm suspension to a volume fraction of 0.08%. The 
non-index-matched beads were used as probes of the thermal fluctuations in our samples. The 
suspension was then injected in a square chamber (3x3 mm 2 ), made of a microscope plate and a 
cover-slip separated by a thin adhesive spacer (250 [im thickness) . The chamber was then sealed 
with glue in order to avoid contamination and evaporation. 

The microgcl particles suspension at T = 30 °C, was found to behave as an equilibrium 'liquid' 
state, characterized by the nearly linear dependency of the probes mean-squared displacement 
with the lag time (see section GO}- It was thus taken as our reference liquid equilibrium state at 
volume fraction denoted $o- 

An estimate of the lower bound for the sedimentation time of the Latex probes over 1 microm- 
eter is 200 s, in the reference liquid of viscosity equal to 6 mPa.s (taking the water density as 
the suspension density). The sedimentation is thus slow enough to allow us recording during few 
minutes. Generally, the sample is left on the microscope at T = 30 C during two minutes for 
recording and then quenched at low temperature to reach high volume fraction. At high volume 
fraction, in the glass state, the sedimentation of the probe particles is negligible due to the high 
viscosity of the material, and recording for hours is possible. 

In soft spheres suspensions, the volume of the microgcl particles might not be fixed at different 
volume fractions, especially at high volume fractions, and the particle number density is no longer 
equivalent to the volume fraction, the control parameter of colloidal suspensions. Moreover, the 
shape of the deformable microgcl particles is not fixed, and high volume fraction suspensions can 
still be fluid above random close packing since particle deformations allow for rearrangements^. 
Here determining the volume fraction $ in the reference state, with the Einstein's viscosity 
relation, gives us a value larger than unity of poor relevance, that indicates an already compressed 
state, and extra shrinkage of the particles. Standard methods, such as counting (nearly index- 
matched particles with the solvent) or weighting, were not possible. In the following, because we 
are interested in the phase behaviour with increasing or decreasing the control parameter, the 
volume fraction of the suspensions at equilibrium at T = 30 ° C will be expressed in terms of 
<f>o at T = 30 ° C. Increasing volume fraction will be achieved by quenching the sample from the 
reference state at <&o and T = 30° C, resulting in a new couple ($, T) and the corresponding 
phase behaviour will be explicited. 

2.2 Dynamic Light Scattering 

Dynamic Light Scattering (DLS) experiments were performed on low concentration suspensions of 
microgcl particles. The set-up is composed of a MELES GRIOT Hc-Ne laser source, operating at 
the wavelength A = 632.8 nm, a Brookhaven Instrument Corp. automatic goniometer furnished 
with a BI-9000AT digital autocorrelator. 

In our slightly polydisperse suspensions, analyzing the time dependence of scattered light was 
shown to be less susceptible to experimental fluctuations than analyzing the angular dependence of 
the scattered light^I. The recorded autocorrelation functions was analysed using a William- Watt 
fit of the first scattering function : 

9i(t)=e-^ 

where r is the decay time and j3 takes into account the slight polydispersity of the sample. The 
exponent j3 was found between 0.8 and 1, corresponding to a polydispersity between 5% and 15%. 

Here, DLS experiments were also performed at different scattering angles q for even more 
dependable results. The decay time r was shown to scale with the scattering wavenumber q 
according to - oc q 2 , indicating a diffusive behaviour of the scatterers (Fig. [T] inset). For each 
temperature, the mean hydrodynamic diameter of the particles Dh could then be deduced from 
the relation : 

1 fc B T 2 

t 37rr?(T)Djj 

where kg is the Boltzmann constant, T the temperature, and rj(T) the solvent viscosity. Figure [1] 
shows the monotoneous decrease of the mean hydrodynamic diameter as a function of temperature. 
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Figure 1: Mean hydrodynamic diameter of the pNIPAm particles as a function of temperature, 
measured with dynamic light scattering experiments. The polydispersity was found to vary be- 
tween 5% to 15%, depending on the realization. A sharp decrease of the microgel particles size 
was observed around T = 32 °C. Their diameter was multiplied by 5 from the fully collapsed state 
at high temperature to the fully swollen state at low temperature, reflecting the small degree of 
cross-linking. Inset: linear relation between the inverse T of the correlation function decay time r 
and q 2 , here at T = 28 °C. The linear dependency validates that only the diffusion of the particles 
is probed, and allows us to determine the particles diameter. 



A sharp decrease of the particle size occurs around T = 32 °C. 

In these DLS experiments, we are aware that we are not in the conventional range of radii 
R for using DLS for particle sizing, since Mie scattering theory applies (qR ~ 1), and internal 
fluctuations of the particles might be probed. However, the good agreement of the 1/r ~ q 2 
law indicates a diffusive behaviour of the microgel particles, giving evidence that the internal 
fluctuations do not contribute to the scattered intensity. We also checked that we could obtain a 
good estimate of the hydrodynamic diameter for calibrated polystyrene microbeads of 1 /im and 
1.8 ^tm in diameter, within a 5% error, in the same range of concentration. We are thus confident 
in our results, within an error bar of 5 to 10%. 

2.3 Particle Tracking 

The thermal motion of the tracers was recorded using a inverted Lcica DM IRB microscope with 
a xlOO oil immersion objective, coupled to a CMOS camera (Phantom v9). The camera was 
typically running at 20 frames per second (fps) in the more concentrated states, and at 500 fps in 
the low viscosity suspensions, with a 56 x 56 /im 2 field of view. Most of the records were 2 • 10 3 
frames long. Depending on the realization, 10 to 70 probe thermal motions could be recorded 
in the same field of view. For reliable analysis of the Brownian motion, particular attention was 
brought to record the motion far from the rigid wall imposed by the glass slides. 

A home-made image analysis software allowed us to track the bead positions x(t) and y(t) close 
to the focus plane of the objective. For each probe i, the time-averaged mean-squared displacement 
(MSD), (Arf (t))? = {[xi(t' + t)- x^t')] 2 + [y^f + t) - yi(t')] 2 ) t , was calculated, improving the 
statistical accuracy. When necessary, an ensemble-averaged MSD over all the probes, denoted as 
(Ar 2 (t)}t> : i, could be computed. In the following, (Ar?(t))f will refer to the temporal-averaged 
MSD for one probe, and (Ar 2 (t))t>,i to the ensemble and temporal averaged MSD. 
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In our aging experiments, the recording time was short compared to the aging dynamics, which 
means that no significative aging was observed during the recording. The system could then be 
considered to be in a quasi stationary state, at a given aging time, during the recording. The 
resolution limit of our set-up was tested using a highly concentrated (95% wt) glycerol solution, 
and was found to be 15 nm. 

Probability Distribution Functions (PDF) of the tracer displacements (self part of the van 
Hove correlation function) were also computed for each tracer : pi (Ax,i) — (6 (Ax — Axi (t)))t>, 
and for the overall sample : p(Ax,t) = (S (Ax — Axi (t)))t',i, for different lag times t, where Aa; 
denotes the elementary one-dimensional displacement along the x axis of the field of view and S 
refers to the delta function. In this study, the system is isotropic, and the behaviour was found 
to be quantitatively the same on both x and y axis. 

2.4 Quenching procedure 

One important consideration for the investigation of the aging dynamics in glasses is the prepa- 
ration and reproducibility of the initial glassy state. In colloidal glasses, one has to quench the 
system from the liquid state by increasing rapidly the volume fraction. In pNIPAM microgel par- 
ticles suspensions, such a possibility is offered by varying the temperature. This can be achieved 
rapidly because the swelling / deswelling rate is of the order of the millisecond—. After a quench, 
microgel particles which are swollen in a few seconds, generate stresses through the sample, that 
relax with time. 

The microscope objective temperature could be set between 20 °C < T < 35 °C, using a 
Bioptechs objective heater coupled with a cooling device. The sample, kept at high temperature 
on a Peltier plate at T = 30 °C, was suddenly placed on the objective, set at a lower temperature. 
The sample temperature was then controlled through the oil immersion objective in contact. The 
equilibration of the temperature sample with the immersion oil by thermal diffusion was achieved 
in a few secondf0- 

All the quenches described here were performed from our equilibrium 'liquid' state at $ = 4>o, 
T = 30 °C, resulting in a new couple describing the sample ($, T). The quench defines the origin of 
aging times t w . At t w = 0, the volume fraction of the suspension increases in a few milliseconds, 
and for sufficiently high volume fractions reached after the quench, a glass was formed out of 
the initial equilibrium liquid. Aging, due to structural relaxations in the deformable soft spheres 
suspension, was observed. 

The relative size of pNIPAm particles to probe particles changes by 20 %, through the quench. 
As a result, it is not the same in the starting liquid state and in the corresponding glass state. 
Moreover, the temperature change slightly affects the particles softness. As for the polydispersity, 
it does not change much between T = 20 ° C and T = 30 °C. As long as these quantities (relative 
size of pNIPAm particles to probe particles, particles softness, polydispersity) remain constant 
during the aging process, after the quench, it will not affect quantitatively our results. It would 
be an issue to determine quantitatively the phase diagram of such a suspension. 

3 Results and discussion 
3.1 Aging dynamics in the glass 

The aging dynamics of the microgel particles suspension after a quench in temperature, was 
investigated by looking at the Brownian motion of the tracers. The mean-squared displacement 
(MSD) of the probes as a function of the lag time, is displayed in Fig. O in the reference 'liquid' 
state at $ = $ U5 T = 30 °C, and after quenches of different amplitudes. 

We checked that no heat convection could take place in our samples. In this geometry, the Raylcigh number 
for water, the lower limit viscosity we used, is 1 for a reasonable temperature difference of 10 °C, which is three 
orders of magnitude below the threshold of thermal convection. 
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Figure 2: Mean-squared displacement of the tracers (1 /xm in diameter) as a function of the lag time 
in the pNIPAM microgel particles suspension, for quenches of different amplitude. Starting from 
an equilibrium liquid state at $ = $o, T = 30 °C, the sample is quenched down to T = 25, 22 and 
21 °C, from the left to the right, (a): at T = 25 °C, the suspension exhibits an equilibrium visco- 
elastic behaviour. No aging was observed over more than one hour. The straight line represents 
the slope 1 corresponding to a purely diffusive behaviour, (b): at T = 22 °C, the sample enters 
a glassy state and slowly ages over several hours, (c): for a deeper quench at T = 21 °C, the 
aging dynamics was found to be faster. (b,c): upon aging, the MSD decreases continuously, and 
develops a relatively smooth 'plateau' for large lag times t, while it is still diffusive for short lag 
times. Moreover, the deeper the quench in the glassy state, the faster the glass aging dynamics, 
i.e. the faster and the deeper the plateau develops. 

In the reference equilibrium state, the MSD of the probes is nearly diffusive, characterized by a 
nearly linear dependency with the lag time, and no aging could be observed over more than 4 hours. 
All the quenches were performed from this 'liquid' state at $ = $o, T = 30 °C. The suspension 
viscosity was estimated with the Stokes relation r\ = AksTt / (67ri?(Ar 2 (t))f i) = 6.1 ± 1.1 mPa-s, 
where R is the tracer radius and fcs the Boltzmann constant. 

For a quench from T = 30 °C to T = 25 °C, the suspension, still at equilibrium, becomes 
slightly visco-elastic, characterized by a non linear increase of the MSD with the lag time; no 
aging was observed over more than one hour (Fig. EJa)). For larger amplitudes quenches - down 
to T = 22 °C and T = 21 °C - the suspension was found to fall out-of-equilibrium and to age over 
hours, indicating that the sample has entered a glassy state, as shown in Fig. [U(b,c). Upon aging, 
the MSD decreases continuously, and develops a subdiffusive 'plateau' at large lag times t, while 
it is still diffusive for short lag times. Moreover, with increasing the quench amplitude, the aging 
dynamics is faster and the suspension goes in a deeper glass state. 

The underlying corresponding picture for the MSD data in the glass state, is the following : the 
first diffusive regime, at short times, is associated with the motion of the tracer in the cage formed 
by the surrounding microgel particles. The second regime - smooth plateau becoming deeper 
upon aging - reflects the elastic response of the cage at the probe length scale. Upon aging, the 
cage becomes stiffer due to structural relaxations of the soft deformable spheres, and the tracer 
gets more trapped. Since the pNIPAm particles are soft, the probe can penetrate the shell of the 
surrounding microgel particles, and the plateau is a 'smooth' one. In these experiments, the cross 
over to a diffusive motion of the tracers at larger lag times (cage rearrangement) was not observed, 
indicating that the tracers remain trapped on the observation time scale. 

Here, the glass transition was dchncd with respect to the out-of-equilibrium regime experienced 
by the suspension above a certain threshold volume fraction. This is reasonable since no crystal 
is expected because of the suspension polydispcrsity. The onset temperature below which the 
suspension enters a glassy out-of-equilibrium state, was found between T = 24 °C and T = 22 °C 
(corresponding to 1.5 $0 < $ < 1.6 $o)- Given the relatively small dependence of the micro- 
gel particles size as a function of temperature in this range (see Fig. Q]), the 'glass transition' 
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Figure 3: Characteristic times r versus aging time t w after quenches at T = 22 °C and T = 21 °C 
, recorded over 5 hours. Here, the characteristic time was defined as the lag time needed to reach 
a certain value I 2 of the MSD. The characteristic time was found to increase exponentially with 
aging time, such as r ~ exp(kt w ). 



temperature could not be determined with more accuracy. 

From the MSD data, a characteristic time r associated with the probes diffusion was extracted. 
This characteristic time was defined as the lag time necessary for a probe to diffuse over a length- 
scale I, such as (Ar 2 (t))t>,i — I 2 - Fig. [3] shows the evolution of the characteristic time r with 
aging time, for two different values of I 2 - 2.5 10 _3 /xm 2 and 5.0 10 _3 /im 2 - after a quench from 
T = 30 °C to T = 22 °C. On these length scales, typically of the order of one tenth of the probe 
and pNIPAm particle radius, the characteristic time grows exponentially with aging time by 2 
decades as t ~ exp(kt w ). Another definition involving Q(1,t) = (e~ Ar ;( r )/' y i t , = 0.6, where 
Q(1,t) is a correlation function equivalent to the intermediate scattering function, gave totally 
similar results^. 

Such an exponential increase of the relaxation time was already observed in other soft glasses 
and gels, often coupled to a transition to an asymptotic aging regime (r ~ t w )— ~— . Here the 
characteristic time we define is not a relaxation time in the usual definition. The asymptotic aging 
regime might develop at later aging times t w , on the considered length scales. Experimentally, no 
cross over from the exponential growth to the full aging regime could be observed, by changing the 
quench depth to accelerate the aging process or measuring the characteristic time over 7 hours. 

In these experiments, we probe the dynamical behaviour of the microscopic environment at the 
scale of the pNIPAm particles, that might be very different from larger scale response. The radius 
of the probes being of the order of the radius of the pNIPAm particles, there is no scale separation 
between the probes and the environment constituting elements. We believe that the probes reflect 
the dynamics of the surrounding pNIPAm colloidal particles, even though the relation between the 
probes and microgels dynamics is not straightforward and depends on a number of parameters, 
especially the chemistry of the probes, and the probe size^ 3 -. 
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Figure 4: Rescaled Probability Distribution Functions (PDF) of the probes displacement as a 
function of the rescaled displacement, (a) in the diluted suspension at equilibrium at $ = $o, 
T = 30 °C, for a lag time t = 0.1 s. In abscisse, the unity corresponds to Ax = 0.6 R, where R 
is the probe radius; (b) in the suspension at T = 30 °C, for a lag time t = 0.15 s. In abscisse, 
the unity corresponds to Ax = 0.38 R; (c) in the glass state after a quench from T = 30°C to 
T = 22 °C, at aging time t w = 60 min, for a lag time t = 1 s. In abscisse, the unity corresponds to 
Ax = 0.42 R. The black curves are Gaussian fits Y = exp (-X 2 ) I spH. In the diluted suspension, 
the PDF is a Gaussian, as expected in a purely viscous liquid. In the suspension at T = 30 °C, 
the PDF slightly departs from a Gaussian at large displacements. In the glass state, the PDF 
exhibits a central Gaussian part for short displacements, and an exponential like decay for large 
displacements, observed at any lag time and aging time. 

3.2 Spatially heterogeneous dynamics 

3.2.1 The spatial distribution of local Gaussian dynamics behind non Gaussian van 
Hove functions 

We now analyze the probe thermal fluctuations by looking at the Probability Distribution Func- 
tions (PDF) of the probes' displacements (or van Hove functions). Near the glass transition, they 
have been shown to exhibit non Gaussian shapes, usually assigned to 'dynamic heterogeneity'. 
Here, we analyze the ensemble- averaged and individual PDF of the tracers. 

Figure 0] shows rescaled ensemble-averaged PDF *J (Ar 2 (t)) t '.i p(Ax, t) as a function of the 
normalized displacement Ax/y/ (Ar 2 (i))t'.i in the equilibrium liquid state at different volume 
fractions, and in the glass state. Experiments performed in the diluted suspension at <I> = 0.1 <&o, 
T = 30 °C (dilution factor 10, purely viscous, viscosity 1.85 ± 0.15 mPa.s), shows a Gaussian 
PDF, as expected in a viscous homogeneous fluid (Fig. QJa)). In the equilibrium reference state 
at $ = <f> 5 T = 30 °C, the ensemble- averaged PDF were found to be nearly Gaussian at any lag 
time (Fig. H[b)), while in the glassy state, they are Gaussian for short displacements, and exhibit 
an exponential like decay for large displacements, at any lag time allowing for a good statistics 
(t < 2 s), and aging time (Fig. BJc)). As shown here, the shape of the ensemble-averaged PDF 
changes with increasing volume fraction : it goes from a pure Gaussian at low volume fraction to 
a PDF with a Gaussian central part for short displacements, and a broad exponential like decay 
at large displacements, at higher volume fraction. Moreover, the tail of the distribution becomes 
more pronounced when approaching the glass transition (increasing volume fraction). 

We now consider the individual motions of the probes through the sample. In the glass state, 
where the ensemble-averaged PDF were found to depart from a Gaussian, the spatial distribution 
of the MSD is wide (Fig. Ufa)). The MSD of the individual probes are drastically different from one 
probe to another, even at shortest lag times where the statistics is the highest. The corresponding 
individual PDF of the probes' displacements are Gaussian, pi(Ax,t) ~ exp (— (Ax(t)) 2 / Ai), at 
any lag time allowing for a good statistics (t < 2 s) and aging times t w . Their variance Aj were 
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Figure 5: (a): MSD as a function of the lag time for different probes (open symbols) and the 
ensemble-averaged MSD over the probes (large filled symbols), in the glass state at aging time 
t w = 60 min, after a quench from T = 30°CtoT = 22°C.A broad spatial distribution of MSD was 
observed for each t. In inset, the corresponding ensemble-averaged PDF is plotted, (b): Individual 
PDF of the displacements for three probes in the glass state after a quench from T = 30 °C to 
T = 22 °C, for a lag time t = 0.05 s, at aging time t w = 60 min. For each probe, the individual 
PDF is a Gaussian with a different width, resulting from a different MSD. (c): the narrow spatial 
distribution in the diluted liquid state at $ = 0.1 $0, T = 30 °C, with the corresponding Gaussian 
ensemble-averaged PDF in inset, (d): the individual particles' displacements distributions are 
Gaussian, with the same widths. 
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found to be widely distributed through the glass sample (Fig. Efb)). The resulting ensemble- 
averaged PDF over all the probes, in the top inset and already shown in Fig. He), exhibits a 
central Gaussian part, and an exponential like decay at large displacements. By comparison, in 
the diluted equilibrium liquid state (<£> = 0.1 $o, T = 30 °C), the MSD are narrow distributed 
(Fig. [He)). The individual PDF of the probes are Gaussian and similar in width (Fig. [5{d)), and 
the resulting ensemble-averaged PDF, shown in the bottom inset and in Fig. SJa), is a Gaussian. 
In the reference equilibrium state (not shown here, $ = $07 T = 30 °C), the distribution of the 
individual MSD through the sample is clearly narrower than in the glass state, but wider than 
in the diluted liquid state. The individual PDF are Gaussian with slightly different widths, and 
the corresponding PDF slightly departs from a Gaussian at large displacements, as shown in Fig. 
lib). 

In the Appendix, we analyze the ensemble-averaged and individual Cumulative Distribution 
Function (CDF) of the probes' displacements in the liquid state and in the glass state, in order 
to clearly identify the shape of the PDF at short and large displacements. Using the CDF to 
graph the displacements distribution of sample data and try to identify it with respect to known 
distributions avoids problems of bin size and bin location. We could validate that the ensemble 
averaged PDF is Gaussian for short displacements, and that its decay at large displacements is 
well adjusted by an exponential over at least 3 decades in the glass state, while local PDF remain 
Gaussian (see Appendix). 

In our system, the individual distributions of the probes' displacements were found to be 
Gaussian, whether at equilibrium or out-of-equilibrium, even though non diffusive, in the con- 
sidered range of lag times 0.05 < t < 2 s, at any volume fraction or age. They are given by 
Pi(Ax,t) ~ exp (— (Ax(i)) 2 /Aj), with a variance A, which becomes widely distributed through 
the sample when approaching the glass transition. The variances Aj characterize the local dy- 
namical environment ^. They can also be considered as a measure of the 'probe mobility' within 
the suspension, since the larger the variance Aj, the more 'mobile' the probe is. 

The scenario behind the non Gaussian van Hove functions is a spatially heterogeneous dynamics 
of the suspension. At low volume fraction, the spatial distribution of the probes' mobilities is 
narrow, revealing an homogeneous suspension, and the ensemble-averaged PDF, p(Ax,t), is a 
Gaussian. With increasing volume fraction, the spatial distribution of the mobilities gets wider, 
revealing a spatially heterogenous dynamics. A sum of local Gaussians pi of different widths is no 
longer a Gaussian. As a consequence, the corresponding ensemble-averaged PDF deviates from a 
Gaussian. The observed non Gaussian tails of the ensemble-averaged p(Ax, t) are thus the result 
of the wide spatial distribution of local dynamical environments Aj through the sample. 

Spatial heterogeneities already appear below the glass transition in the equilibrium state at 
$ = $0, T = 30 °C. This suggests that the suspension is already in a 'supercooled' regime, where 
a qualitative change in the character of particles motion occurs. The crowding of the suspension 
is high enough to induce small heterogeneity of dynamics in the sample, but not enough to allow 
aging, out-of-equilibrium fall, or jamming to occur—. 

We now characterize the shape of the distributions f (Aj) of the individual probes' mobilities, 
at equilibrium and out-of-equilibrium. Fig. [5] shows the spatial probability distribution functions 
of the Aj, defined such as (Arf {t)) t > = Aj(t). In the reference liquid state ($ = $0, T = 30 °C) 
and in the glass state, the distribution is a Gaussian for short diffusion coefficients, and exhibits 
an exponential decay for large diffusion coefficients (Fig. EJb,c)). In the diluted liquid state 
(<& = 0.1$o, T = 30 °C), the distribution is centered around a mean diffusion coefficient, and no 
exponential decay was observed for large diffusion coefficients (Fig. [S[a)). Given a good statistics 
for the calculation of Aj(i), the shape and behaviour of the distributions remains similar for various 
lag times t: the distribution -P(Aj) was found to develop an exponential tail when approaching 
the glass state. In these experiments, no clear dependence of the distribution shape was found 
with the aging time. 

It is now straightforward to show that, given an exponential tail for the distribution P(Aj), 



2 Note that the variance Aj of the individual PDF is by definition directly related to the individual MSD, by 
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the resulting ensemble-averaged PDF p (Ax, t) will exhibit an exponential decay at large dis- 
placements. We consider the ensemble-averaged PDF as the sum of local Gaussian dynam- 
ics p(Ax,i) = J dAi P (A,-) cxp (-(Ax(t)) 2 / A;) / \JnAi, with a spatial distribution P (Aj) = 
exp (~^'/ A °) i a t l ar g e A». The larger the Aj, the more it contributes to the integral, so we can 
write : p(Ax,t) ~ / dA, cx p(-^/a ) exp (-(Ax(t)) 2 / A,) /V^A", which yields, using a method of 
steepest descent, a leading term for the ensemble averaged PDF at large displacements p (Ax, t) ~ 
exp (— 2|Aa;(t)|/\/Ao). The exponential decay of the distribution P(Aj), at large Aj, results in 
an exponential decay of the ensemble-averaged PDF. 

The exponential decay of P(Aj) reveals that a small fraction of tracers are intrinsically faster 
than average. Given that all the individual PDF are Gaussian, the fast particles are responsible 
for the large displacements in the ensemble-averaged PDF. Here, the exponential decay of the 
ensemble-averaged PDF is thus the result of the existence of faster mobility zones than average 
on our observation time scales, rather than jumps experienced by the tracers from time to time. 
Moreover, in these experiments, no evidence of temporally heterogeneous dynamics was detected 
on the considered recording times, in the sense that fast regions did not switch from rapid to slow 
dynamics over time (and vice versa). 

We have characterized the statistical distribution of dynamical environments, and show that it 
can explain the observed exponential tail of the van Hove functions, in the concentrated regimes 
(supercooled and glass). We analyze quantitatively below how the variance cr 2 (t) of this distribu- 
tion P{Ai) evolves when approaching the glass transition, and beyond. 



3.2.2 Spatial heterogeneities when approaching the glass transition, and beyond 

A question that naturally arises now to further characterize the spatial heterogeneities concerns 
the way they evolve, when approaching and entering the glass state. We now look at the variance 
u 2 (t) of the spatial distribution of the Aj(i) = (Ar 2 (t)) t > through the sample, given by : 

a 2 ({Ar 2 (t)) v ) (((Ar 2 (t)) t ,-(Ar 2 (t)) t ,,,) 2 ) t 
(Ar 2 (m ti (Ar 2 (i))^ 

where (Ar 2 (t)) t ' refers to the temporal-averaged MSD for one probe, and (Ar 2 (t)) t \i to the 
ensemble and temporal averaged MSD. This variance gives a measure of the intensity of spatial 
heterogeneities in the sample. The larger the variance a 2 , the more distributed the variances 
Aj of the individual tracer motion through the sample, and the more the sample is spatially 
heterogenous. 

Figure [7] shows the variance a 2 of the spatial distribution of the (Ar 2 (t)) t > through the sample, 
as a function of the volume fraction, before and after entering the glass state. Starting from a very 
diluted liquid state, the variance was found to increase with increasing volume fraction. Given 
a good statistics for the calculation of A,(i), the variance <J 2 (t) was found to be independent of 
the lag time t, and Figure [7] unchanged. The water represents the reference where no spatial 
heterogeneities are expected, and sets the sensitivity in our experiments. In thermosensitive 
suspensions, the value for a 2 (t) was found to be larger than the experimental sensitivity, which 
indicates that spatial heterogeneities are present. Figure [JJ shows that the suspension dynamics 
becomes more and more spatially heterogeneous when approaching the glass transition. As seen 
in section 13.2.11 the dynamics is spatially homogeneous at low concentration ($/<&o — 0.1, T = 
30 °C), and becomes spatially heterogeneous, already below the glass transition (Q/$o = 1, T = 
30 °C). 

When entering the glass state (above the vertical line) , the dynamics remains strongly hetero- 
geneous but also widely distributed. The behaviour of spatial heterogeneities in the glass state 
during aging is shown in Fig. [FJ where we plot the variance a 2 as a function of a given MSD, 
(Ar 2 (t)) t ^i, at t = 2 s. We take this MSD as a measure of how deep the suspension is in the glass 
state after a quench, since this allows us to compare quenches with different final temperatures 
and different realizations of these quenches, from the same initial equilibrium state. For a given t 
in the diffusive regime or in the plateau, the variance tends to increase with decreasing amplitude 
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Figure 6: Probability distribution functions of the variances A h defined such as (Arf(t)) t > = A;(f). 

(a) in the diluted liquid state at (T = 30 °C, $ = 0.1$ ), over 1200 probes, lag time t = 6 10" 3 s; 

(b) in the reference liquid state at (T = 30 °C, $ — $0)7 over 1150 probes, lag time t = 3.75 10 -2 
s; (c) in the glass state at T = 22 °C and t w = 45 min, over 1220 probes, lag time t = 6 10~ 2 s. 
(b,c): in the reference liquid state and in the glass state, the distribution is a Gaussian for short 
diffusion coefficients, and exhibits an exponential decay for large diffusion coefficients. No clear 
dependence of the distribution shape was found with the aging time, (a): in the diluted liquid 
state, the distribution is centered around a mean diffusion coefficient, and no exponential decay 
was observed. 
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of the (Ar 2 (t)) t \i- This suggests that the deeper the suspension in an arrested state, the more 
spatially heterogeneous the dynamics. 

This work establishes the link between the non Gaussian, exponentially tailed van Hove func- 
tions and a spatial distribution of homogeneous dynamic environments in a soft colloidal glass. 
The intensity of the dynamic heterogeneities is characterized when approaching the glass transi- 
tion. In aging experiments, reproducibility in the out-of-cquilibrium dynamics after a change in 
the external control parameter, is a major experimental issue. In our experiments, we noticed that 
taking the aging time t w as a parameter of the 'glass state' gives qualitatively contradictory results 
for a 2 (t w ), from one realization to the other. However, choosing a given MSD, (Ar 2 (to))t>,i, as 
a measure of the 'glass state' leads to the same qualitative behaviours for a 2 as the suspension 
goes deeper in the glass state. The dispersion of the data in Fig. [5] certainly comes from the poor 
statistics consecutive to the use of tracers instead of the microgcls themselves. More systematic 
experiments in the aging regime, improving the statistics, in suspensions with particles of various 
softness, are currently under investigation. 

Dynamic heterogeneities in the aging regime have been investigated in a few studie a 22 ' 23 i 44 , 
although results in this regime might help to discriminate between theories^. Since the particles 
arc almost index matched and impossible to visualize unambiguously on the microscope, the 
characterization of the shape and size of the regions of homogeneous dynamic environments were 
beyond the scope of this study, as well as an attempt to question a link between the structure and 
the dynamic heterogeneities^. Here, the dynamic heterogeneity is seen as a spatial heterogeneity 
on the recording time scale, which was chosen to ensure quasi static conditions. Complementary 
work will investigate the temporal fluctuations of the dynamics and the relation between aging 
and temporally heterogeneous dynamics. 



4 Conclusion 

We have investigated the aging dynamics and dynamic heterogeneity of a soft thermosensitive 
microgels suspension, by looking at the thermal fluctuations of tracers. By varying temperature, 
the volume fraction could be tuned, and the characteristics of the equilibrium liquid suspension 
and the out-of-cquilibrium glass state during aging, investigated on the same sample. Aging was 
observed over several hours by quenching the suspension from an initial liquid state at high tem- 
perature. With increasing the quench amplitude, the aging dynamics was found to be faster, and 
the system to enter a deeper glass state. Spatially heterogeneous dynamics was observed to de- 
velop below the glass transition, and to persist in the glass state. Individually, the probes always 
exhibit a Gaussian PDF, whether at equilibrium or in the glass state, revealing a locally homoge- 
neous dynamical environment. With increasing volume fraction, the local dynamical environment 
was found to change significantly from one probe to another, sign of a spatially heterogeneous 
dynamics. In the supercooled regime, characterized by a non aging heterogeneous system, and in 
the glass states, the statistical distribution of these dynamical environments was shown to decay 
exponentially, indicating that the dynamics is dominated by a small fraction of faster mobility 
zones than average, rather than jumps experienced by the tracers from time to time, on the ex- 
perimental timescales. We show that these distributions with exponential tails can explain the 
observed exponential tails of the van Hove functions. Finally, we analyze the variance of these 
distributions, as a function of the volume fraction, and as a function of the age in the glass state. 
It gives a measure of the spatial heterogeneity intensity in our samples. With increasing volume 
fraction, spatially heterogeneous dynamics was found to amplify. In the glass state, it tends to 
increase as the suspension gets more arrested. 
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Figure 7: Variance of the spatial distribution of the individual MSD as a function of the relative 
volume fraction, at equilibrium and in the glass states. The relative volume fraction is defined with 
respect to the volume fraction $o of the reference equilibrium state at T = 30 °C. Given a good 
statistics for the calculation of Aj(t), the variance a 2 {t) was found to be independent of the lag 
time. The water represents the reference where no spatial heterogeneities are expected, and sets the 
sensitivity in our experiments. In thermosensitive suspensions, the value for u 2 (t) is larger than the 
sensitivity, which indicates that spatial heterogeneities are present. At equilibrium, the variance 
increases with increasing volume fraction (roughly exponentially for large volume fractions, as a 
power law for low volume fractions, see inset), indicating the suspension becomes more and more 
spatially heterogeneous when approaching the glass transition. In the out-of equilibrium states 
(above the dashed vertical line), the suspension remains strongly spatially heterogeneous and also 
widely distributed (the different sets of points correspond to different quenches; within a quench, 
the different points correspond to different aging times), (inset): corresponding log- log plot. 
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Figure 8: Variance of the distribution of the individual MSD in the glass state, after a quench 
from T = 30 °C to T = 22 °C, as a function of a given MSD, at a lag time t = 2s. At T = 30 °C, 
the suspension viscosity iss 6.0 ± 0.3 mPa.s. Although widely distributed, the variance a 2 tends 
to decrease with the mean MSD, suggesting that the more the suspension is arrested the more the 
dynamics is spatially heterogeneous. The dashed line is a guide to the eye. 



5 Appendix 

5.1 Cumulative Distribution Function 

The Cumulative Distribution Function (CDF) C (x) of a random variable (e.g. the displacement 
of a probe particle during a fixed lag time, Ax) is the probability that this random variable takes 
a value smaller or equal to x : C (x) = p (Ax < x). The CDF of a random variable is deduced 
from its PDF p by integration : C (x) — J_ dx'p (x 1 ). A schematic representation of a typical 
CDF is presented below : 



0,5 



0, 




It monotonically increases from to 1, is antisymmetric around the variable median (0 in our 
case) and has a typical width of the order of its standard deviation (square root of the MSD in 
our case). 

Our interest here is to know the detailed shape of the tails of the CDF, for large positive and 
negative displacements. Since the situation is symmetrical, we focused on positive displacements 
only. In order to zoom on what happens as C (x) reaches 1, the tail is stretched by computing, 
for x > 0, the quantity : 

y = -ln(2(l-C(x))), (1) 
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which diverges as C (x) goes to 1 (x to +00). For the negative tail, we computed y 1 = ln(2C (x)) 
and found results in quantitative agreement with those of the positive tail. 

Two particular cases were of great interest to us : the case of a Gaussian PDF of the displace- 
ments and the case of a PDF of displacements decreasing exponentially for large displacements. 
In the case of a Gaussian pdf of the displacements p (Ax) = exp (— (Ax) 2 /(Ar 2 )) / yj (Ar 2 ), the 
CDF is given by : 

1 + erf (s/v^A^y) 

C gauss (x) = - . (2) 

For large displacements (x >> \J (Ar 2 )), the stretched CDF is given by inserting Eq. [2] in EqQ] 
and developping for large x and yields : 

V gauss = U 2 + - In 7T + 111 U + O (u~ 2 ) , (3) 

where u — xj y (Ar 2 ). Thus a quadratic variation of y with x at leading order is reminiscent of a 
Gaussian shape of the PDF of displacements. We fitted variable y in the Gaussian cases with its 
full expression (Eq. [5] inserted in EqfT]) rather than the large displacement development. 

In the case of a PDF with exponential tails, if the PDF of displacements follows the equation 
p (Ax) = exp (— |Ax|/A) /A for large Ax, then The CDF of displacements is given for large x by : 

C exp (x) = l-Be-%, (4) 

where B depends on the detailed shape of the PDF at small displacements. The stretched CDF 
is then given for large displacements by : 

y eap = ^-ln(2S). (5) 

A linear divergence of the stretched CDF is thus equivalent to an exponential decay of the PDF 
of displacements. 

Finally, if the stretched CDF is fitted by y = au + b (u = xj (Ar 2 )), it is easily shown that the 
PDF of displacements is given by | exp (—6) exp (—au). If the fit is valid on a range u <G [it ; U\], 
then the PDF decrease exponentially on a number n of decades given by n = log 10 (p(uq) / p(u\)) = 
(a{m -u ))log 10 (e). 

Figure[9]shows the individual and ensemble-averaged Cumulative Distribution Functions (CDF) 
of the probes' displacements in the low concentration liquid equilibrium state and in the glass state. 
In the low concentration liquid state, the individual CDF are narrow distributed, indicating a spa- 
tially homogeneous suspension. All the CDF in the liquid state, individual and ensemble-averaged, 
could be adjusted by a parabolic curve indicating a Gaussian behaviour of the corresponding PDF. 
In the glass state, the individual CDF are widely distributed and still described by a parabolic 
function, as seen in Figure fright. Their ensemble-average results in a linear function, giving 
evidence of an exponential decay of the associated PDF over 3 decades. 

In the glass phase at T = 22°C, we found a = 2.40 ± 0.05 and b = -0.40 ± 0.05, valid over a 
range u <G [0.7; 3.5]. This is equivalent to an exponential decrease of the PDF over 3 decades. 
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